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The density matrix of a two-level system (spin, atom) is usually determined by measuring the 
three non-commuting components of the Pauli vector. This density matrix can also be obtained 
via the measurement data of two commuting variables, using a single apparatus. This is done 
by coupling the two-level system to a mode of radiation field, where the atom-field interaction is 
described with the Jaynes-Cummings model. The mode starts its evolution from a known coherent 
state. The unknown initial state of the atom is found by measuring two commuting observables: the 
population difference of the atom and the photon number of the field. We discuss the advantages 
' of this setup and its possible applications. 
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INTRODUCTION 



Determining the unknown state of a quantum system is the basic inverse problem of quantum mechanics. Given 
the state, one can calculate the expectation value of any observable of the system However, the inverse problem of 
determining the state by performing different measurements is non-trivial. This problem was discussed by Pauli 0] in 
Q-T 1933. His question was whether one can reconstruct the unknown wave-function of an ensemble of identical spinless 
particles via the corresponding position and momentum probability densities. The interest to the state determination 
problem grew considerably since then, and now this is a well-recognized subject 0]. 

The notion of state refers to an ensemble of identically prepared systems and is represented by a Hermitian operator 
with non-negative eigenvalues that sum to one (density matrix) [l| . Thus the N x N density matrix of an ensemble of 
TV-level systems is specified by N 2 — 1 independent real parameters. Since any observable O of the system generates at 
most N — 1 independent probabilities, at least (N 2 — 1)/(N— 1) = N+l measurements of non-commuting observables 
are needed to obtain the unknown state Q. 

Procedures of reconstructing the quantum state from measurements are known as quantum state tomography. 
Recently they found applications in quantum information processing [j| . For example, in quantum cryptography one 
needs a complete specification of the qubit state @ . For the simplest example of a spin-i system the state is described 
by a 2x 2 matrix. According to the above argument, one has to perform 3 incompatible measurements for the unknown 
state determination, e.g., measuring the spin components along the x-, y- and z- axes via the Stern-Gerlach setup. 
However, during the measuremental procedure of each component one looses the information about the two other 
r^"* , components, since the spin operators in different directions do not commute. Thus, to determine the state of a spin-i 
system, one needs to use repeatedly three Stern-Gerlach measurements performed along orthogonal directions. 

However, the state can be characterized indirectly via a single set of measurements performed simultaneously on 
the system of interest and an auxiliary system (assistant), which starts its evolution from a known state [§]. This 
can be done by letting them interact for a specific time and measuring commuting observables of the overall system 
(system+assistant). In particular, it has been shown that one can determine the unknown state of a spin-i system 
with a single apparatus by using another spin-^ assistant The idea was recently implemented by Peng et. at Q 
who used pulses to induce the proper dynamics of the interaction between a spin-i system and its assistant. They 
verified the initial state of the system obtained from this procedure with the results of the direct measurement of the 
three components of the spin vector of the system. 

Our present purpose is to determine the unknown density matrix of an ensemble of two-level systems (atom or 
spin) via interaction with a single mode of the electromagnetic field. The atom-field interaction is studied within 
the Jaynes-Cummings model [T(| • We show that the unknown state of the spin can be completely characterized by 
measuring two commuting variables: the population difference of the atoms a z and the photon number of the field 
a^a. This measurement supplies three averages: (o~ z ), (aJa) and (o~ z aJa), which will be linearly related to the elements 
of the initial density matrix of the ensemble of the two- level atoms. (Note that since a z and a) a commute, (a z a^a) 
is recovered from the measurement data of o~ z and a' a via the number of coincidences.) 

In section [5] we will give a brief introduction on the Jaynes-Cummings model and its properties. In section [3] the 
model is used to determine the state of the two-level system. Section 13. Al discusses imperfection of the proposed 
scheme due to the noise in selecting the measurement time. In section 0] we discuss how to reconstruct the unknown 
density matrix approximately given an incomplete measurement data. The solution of this task amounts to a di- 



2 



rect application of the classical Maximum Likelihood setup, because our scheme operates with the measurement of 
commuting observables. We conclude in section [5j 

2. THE JAYNES— CUMMINGS MODEL 

The Jaynes-Cummings model ( JCM) [l(| has a special place in quantum optics and atomic physics [ill, E2, EH • 
This model describes the interaction of a two-level atom with a single mode of electromagnetic field, and it was 
employed by Jaynes and Cummings for studying the quantum features of spontaneous emission. Later on, the model 
generated several non-trivial theoretical predictions — such as collapes and revivals of the ato mic p opulation that are 
related to the discreteness of photons [14, 1 1 51] — which were successfully tested in experiments [161 ] . In particular, the 
model explains experimental results on one-atom masers [lH , and on the passage of (Rydberg) atoms through cavities 
[]~7L El, Ell H3] • The JCM is also used for describing quantum correlation and formation of macroscopic quantum 
states. It was recently employed in quantum information theory [U, More recently, the JCM found applications 
in semiconductors [251 ] and in Josephson junctions [24 ] . 

The JCM is in fact a family of models, since the original model of Jaynes and Cummings was generalized several times 
for a more adequate description of the atom-field interaction (e.g., multi-mode fields, multi-level atoms, damping) 
[23L 1261 I27I ]. We shall however study the simplest original realization of the JCM that involves a two- level atom 
interacting with a single mode of electromagnetic field. In particular, we neglect the effects of noise and dissipation. 
This situation has direct experimental realizations [lrl |28| . . For instance with the superconducting microcavities 
one can achieve ~ 0.1s for the average lifetime of the cavity photon, much larger than the typical field-atom interaction 
time ~ 100 - 500^s Hj. 

The quantized mode of the radiation field is described via bosonic creation a' and annihilation a operators. The 
two-level system is mathematically identical to a spin-i, and can be described with the help of the Pauli matrices a x , 
a y and a z . Within the dipole approximation one has the following Hamiltonian for the atom and the cavity mode: 

H = hwa z + hva}a + hg(a' + a)a x (2-1) 

where u> and v are the atom frequency and the mode frequency, respectively, and g is the atom-field coupling constant. 
In the dipole approximation it is given as g = d^/hVeo) 1 ^ 2 , where d is the atomic dipole matrix element, and V is 
the cavity volume. 

We note that in quantum optical realizations of the JCM the coupling constant g is normally much smaller than uj 
and v, e.g., it is typical to have v ~ ui x 10 5 (oj ~ v = 10GHz, g ~ A =10-100KHz). Thus the subsequent reasoning 
based on the interaction representation is legitimate. To obtain from (|2.ip the JC Hamiltonian note that in the 
interaction representation the coupling term reads: 

hg{a) e- wt + ae wt )(<7_ e - *"* + fr+e*"*), (2.2) 

where we introduced raising 

& + = a x + ia y , (2.3) 

and lowering 

0- = a x - ia y , (2.4) 

spin operators. Recall that they satisfy the following commutation rules 

[ff±,^]=Tff±, [&+,&-) =2a g , a+a- +<j_a + = 1. (2.5) 

One now applies to (|2.2p the rotating wave approximation: the atom and field frequencies are assumed to be close 
to each other, and then the factors oc e ±it ^ +w ) in (|2.2[) oscillate in time stronger than oc e ± lt i v -^ . Thus the factors 
oc e ±lt{v+ ^ are ne glected within this approximation and one arrives at the JC Hamiltonian: 

H = hiu<T z + hva) a + hg(a+a + <7_aJ). (2-6) 

We shall denote 



A = LO — v 1 



(2.7) 
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for the detuning parameter. For our future purposes we consider A as a tunable parameter. Within the atom-cavity 
realizations of the JCM, the detuning A can be controlled via the mode frequency v, that is, by changing the shape 
of the cavity. Alternatively, A can be changed via the atom frequency u> by applying an electric field across the cavity 
[30| . Then ui is modified due to the Stark effect. 
The above standard derivation of (|2.6[) is based on small detuning A and weak atom-mode coupling g: 

(2.8) 



A <C min(o;, u), g <C min(w, v). 



Both these conditions are usually satisfied for quantum optical realizations of the JCM. 

There are however situations — especially for the solid state physics applications of the Hamiltonian (|2.1[) — where 
the atom-field interaction constant g is not small. To this end it is useful to know that sometimes the counter-rotating 
terms cx e ± lt (' y + u ') vanish due to specific selection rules [31(, and then the JCM applies in the strong-coupling situation 
as well. 

The Hamiltonian (|2.6p is exactly solvable and the corresponding unitary evolution operator reads [ll[ 



U(t) = coa[ty/<p + g 2 
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where the operator ip is defined as 



The unitarity of U is guaranteed by the identities 



-)(" 



(2.9) 



(2.10) 



sm 



ty/ip + g* 



sin [t y/tp] 



a cos t\J p 



(2-11) 



3. DETERMINATION OF THE ATOM INITIAL STATE. 



Let initially the atom be described by some general mixed density matrix ps'. 

Ps = ^(1 + (v x ) a crx + (v y ) Vy + (vz) a crz), (&i) = tr(CTj ps), i = x,y,z, 

where (<Ti) are the three unknown coefficients of the initial atom state. 

We shall assume that the mode starts its evolution from a coherent state with a known parameter a: 

\a)=e-^/*jr i ^=\n), 



where \a) is the eigenvector of the annihilation operator a, 

a\a) = a | a), 

and where \n) is the eigenvector of the photon number operator a^a, 

a) a\n) = n\n). 



(3.1) 



(3.2) 

(3.3) 
(3.4) 



4 




200 

Interaction time [fis\ 



200 300 
Interaction time [fis\ 



(a)n = 2, A = lOKHz 



(b)n = 2, A = lOOKHz 




2111) 

Interaction time [^s] 



(c)fi = 5, A = lOKHz 



100 200 

Interaction time [jjsj 

(d)n = 5, A = lOOKHz 






ioo 200 : 

Interaction time [^s] 

(e)n = 10, A = lOKHz 



100 200 31 

Interaction time [jjsj 

(f )n = 10, A = lOOKHz 



FIG. 1: The determinant as a function of time for different values of the mean photon number, n = 2, n — 5, and n — 10 and 
detuning parameter, A = 10 and A = 100 KHz. In all cases D has larger maximum value when A = lOOKHz. The coupling 
constant is fixed at g = 50KHz. 



For the average number of photons we have: 



n = (a a' a a) 



(3.5) 



The assumption (13. 2ft on the initial state of the field is natural since these are the kinds of fields produced by 
classical currents [32j |. and also, to a good approximation, by sufficiently intense laser fields. 

Since initially the system and the assistant do not interact, the overall initial density matrix is factorized 

p(0)=ps®\a)(a\. (3.6) 
With the help of the unitary operator (|2.9|) one can calculate the overall density matrix p{t) at time t: 

p(t) = U(t)p(0)U\t). (3.7) 
Then the expectation value of any observable O of the overall system at time t is 

(6) = Tr (p(t)6) . (3.8) 
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In the Appendix A we work out equations for the atom population difference (cx z )t, the average number of photons 
(aJa)t and the correlator of these two observables (a z a'a)t] see (|A1[) - (IA7|) . Expectedly, these three quantities are 
linearly related to the three unknowns parameters (<7 x )o, (&y)o an d (&z)q of the initial atom density matrix: 

B = . (3.9) 




The elements of the matrix M. and of the vector B are read off from (|A1|) - (|A7[) ; see as well Appendix A. They 
depend on the parameter a of the initial assistant state, on the parameters A and g of the JC Hamiltonian, and on the 
interaction time t. Thus, if the determinant of M is not zero, one can invert M and express the unknown parameters 
of the initial atom density matrix via known quantities. Although the elements of M. are rather complicated, the 
determinant itself is much simpler. It takes the explicit form 

00 00 |_,|2(n+m+l) 

# = dctM= 4Ag 2 e- 2 l"' V V 11 , , (n - m) x 



n— m— 



sin 2 (%^) sinf2 m ^ sin 2 (%^) sinf2 n t 



n 2 n n 2 n 



(3.10) 



where O n is the corresponding Rabi frequency 



Q n = V4(« + l)5 2 + A 2 . (3.11) 

Eq. (|3.10jl is our basic result. At the initial time t = 0, the determinant D is zero, since the initial state of the 
overall system is factorized. It is seen in Figs. [TJ(a)-[TJ(f) that for a non-zero detuning, A / 0, the determinant D is 
non-zero for a certain initial period t > 0. (Obviously D is zero when there is no photon in the cavity.) Comparing 
figures Fig. [TJ(a) and Fig. [T](c) we see that although higher initial photon numbers n lead to bigger values for the 
determinant, they cause rapid oscillations in the value of D. This makes the measurement process more difficult. 
(Note in this context that the determinant depends on the absolute value of a.) 

If the average number of photons n = \a\ 2 in the initial state of the field is sufficiently large, D collapses to zero for 
intermediate times; see Figs. [U(e) and[TJ(f). The reason for this collapse is apparent from (|3 . 10[) and has the same 
origin as the collapse of the atomic population difference well known for the JCM [12j. Each term in the RHS of 
(I3.10[) oscillates with a different frequency. With time these oscillations get out of phase and D vanishes (collapses). 
However, since the number of relevant oscillations in D is finite, they partially get in phase for later times producing 
the revival of D, as seen in the Figs.[T](e) and[TJ(f). 

It is seen that D does not depend on separate frequencies u) and v of the two-level system and the field, only their 
difference A = w — v is relevant. This is due to the choice of the measurement basis — see the LHS of (I3.9|) — that 
involves quantities which are constants of motion for g — > 0. Note that D = for A = 0. Thus some non-zero 
detuning is crucial for the present scheme of the state determination. The value of D changes by varying the detuning 
parameter A. Comparing the figures Fig.[T](a) with Fig.[TJ(b), Fig.[T](c) with Fig. ffl(d), and Fig.[TJ(e) with Fig.[T](f) 
one observes that the value of the highest pick of D increases by an order of magnitude when the detuning parameter 
changes from lOKHz to lOOKHz. Note that in Eq. (|3 . 1 0[> for the determinant D the contribution from the diagonal 
n — m matrix elements of the assistant initial state |o;)(a| cancels out. Thus, it is important to have an initial state 
of the assistant with non-zero diagonal elements in the {|n)} basis. 

The basic message of this section is that the determinant D(t) is not zero for a realistic range of the parameters. 
This means that the initial unknown state of the two level system can be determined by specifying the average atom 
population difference (a z )t, the average number of photons (a^a) t , and their correlator (a z a) a)t- In their turn these 
quantities are obtained from measuring two commuting observables: the atom population difference & z and the photon 
number a)a. Having at hand the proper measurement data for these two observables, one can trivially calculate (<r z )t, 
(a'a)t, and find (a z aJa)t via the number of coincidences. 



1. Numerical illustration 



At this point it may be instructive to give two concrete examples on the inversion of the matrix M. in (|3.9p . 
1. Let us ssume that the average number of photons inside the cavity is two n = 2, the coupling constant is g = 50 
KHz, and the detuning parameter A = 10 KHz. Looking at Fig. 1(a) one sees that the determinant is maximal at 
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(approximately) r = 20/is. (Recall that the typical interaction time of a thermal atomic beam with the single mode 
of the field is of the order of 100/ks [l?], H3|-) The elements of the matrtix M. and the vector B are worked out in 
Appendix. Inserting all these numbers into (|Al|) - (|A8jl one obtains 

/ 15.183 5.59578 0.0456968\ 
M- l \&= WKH z = 1-14077 -1.3668 -1.38923 , (3.12) 

V 1 1 j 

/-0.055763l\ 

B\a=wkh z = 2.05576 . (3.13) 
\ 0.0411884 J 

2. For the second example we take a larger detuning: n = 2, g = 50 KHz and A = 100 KHz. The proper 
interaction time t = 300/xs is read off from Fig. 1(b) (interaction time t ss 18 /is gives somewhat smaller determinant; 
see Fig. 1(b)). The numerical calculation of M.~ l and B produces: 

/3.18085 1.23251 -1.15186\ 
M~ 1 \a=wokhz = 5.92052 -4.20194 -4.52702 , (3.14) 

V 1 1 J 

/-0.0707962\ 

B\a=wokh z = 2.0708 . (3.15) 
\ -0.119635 / 

A. Random interaction time. 

We saw above that the success of the presented scheme is to a large extent determined by the ability to select 
properly the interaction time t, since this ultimately should ensure a non-zero (and sufficiently large) determinant D 
(It is clear that a small determinant will amplify numerical errors; see the next section for an example). 

To quantify the robustness of the presented scheme it is reasonable to assume that there is no perfect control in 
choosing the interaction time. To this end let us assume that the interaction time t is a random, Gaussian distributed 
quantity centered at to with a dispersion a. The corresponding probability distribution P{t) of thus reads 

Pit) = — e -(*^o) 2 /(2-). (3 . 16) 
27TCT 

One now averages the determinant D over this distribution, 

oo oo i |2(n+m+l) 

D(t )= 4A g 2 e~ 2 H- V V 11 (n - to) [w(Q n ,n m ;t ) - w(Q m ,n n ;t )} , (3.17) 



n'.ml 

n— m— 



where 



w(fl n ,Q m ;t 
1 



2e-% n ™~ sm[t n m } - e~* sm[t (n m + n„)] - e "f sin[(n m - n n )t ] (3.18) 



It is seen that the oscillations of D(t) turn after averaging to exponential factors e -%-( n m±o„) 2 anc j e -V n m in 
()3.17l |3~T5|) . due to which the averaged determinant D(to) gets suppressed for a sufficiently large "indeterminacy" a. 
This suppression is illustrated in Figs. 2 (a) and 2 (b). 

4. MAXIMUM LIKELIHOOD RECONSTRUCTION OF THE INITIAL SPIN STATE. 

Above we have shown how one determines the initial spin density matrix given the three averages (&z)t, (cfta)t, and 
(<r 2 a)a)t- However, in practice the measurement statistics that leads to the above three averages may be incomplete 
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FIG. 2: The averaged determinant D as a function of to f° r different values of a; see (|3.16fl57I5|l . Here g = 50KHz. 

(due to various noises and experimental imperfections) and we should understand how to reconstruct the initial 
density matrix approximately given the incomplete measurement data. The general question on the approximate 
state reconstruction (given incomplete statistics) is of obvious importance and it got much attention in the standard 
schemes of quantum tomography; see [33[ and references therein. Recall that in these schemes one measures non- 
commuting observables. In this context Refs. [33| propose to generalize suitably the classical Maximum Likelyhood 
(ML) method; see [34j for a detailed discussion on this inference method. In particular, this generalization accounts 
for the fact that the (incomplete) data is obtained from measuring non-commuting observables. 

Once our scheme operates by measuring the commuting variables only, we are going to show that for the approximate 
state reconstruction one now does not need anything beyond the most standard (classical) ML method. Since one 
measures the number of photons and the spin direction along the z-axes (these quantities are represented by the 
operators cva and a z , respectively), the incomplete data in our case means that we are given frequencies v a (m) of 
events, where one registered m photons (to = 0, 1, . . .), and where, simultaneously, the spin component assumed values 
a = ±1. Now recall [34[ that in the ML method the probabilities p a (m) (given the frequencies v a {m)) are obtained 
by maximizing over p a (to) the likelihood function 1 



C\p a (m)]= ^ v a (m) In p a (T 



(4.1) 



a— ±1 m— 



This maximization over p a (m) is to be carried out in the presence of relevant constraints. For our case the initial 
spin density matrix ps must be a positive-definite, normalized matrix; see (|3.6p . Thus we get a single constraint 



'v)o 



(4.2) 



Working out (|3.9[) we write this constraint as a function of the probabilities p a (to) : 

(u - B) T C{u - B) < 1, (4.3) 
where T means transposition, C = (AiAi 1 ")^ 1 , the matrix M. and the vector B are defined in (|3.9|) , and where finally 



VEa=±i E™=o amp a (m)J 
If the constraint (|4.3[) is satisfied automatically, the maximization of C[p a (m)] in (|4.ip produces [34 

Pa(m) = v a (m), 



(4.4) 



(4.5) 



i.e., that the sought probabilities are equal to the frequencies, as one would expect intuitively. However, in general this 
constraint is not satisfied automatically and has to be included explicitly in the maximization of C[p a (m)] over p a (m). 



1 Equivalently, one can mimize over p a (m) the relative entropy 53 a =±i 5Zm=o Ua ( m ) ' n "°^ m ^ 



This measure of distinguishability between 



p a (m) and u a (m) is equal to zero if and only if p a (m) = u a {m) and it has an important information-theoretic meaning [34 



8 



Indeed, looking at (|4.ip and (]4.3[) we may deduce qualitatively that the constraint (|4.3p will be satisfied automatically 
by (|4.5[) . if the frequencies are not very far from the actual probabilities (the ones that would be obtained in the 
perfect experiment) and, simultaneously, the determinant det.A/1 is not very close to zero. 



1. Numerical illustration 



Below we give a concrete numerical example, where the constraint (|4.3[) may or may not be satisfied automatically. 
We take n — 2, g = 50 kHz, A = 100 kHz, and we have chosen the measurement time t = 300/is such that the 
corresponding determinant D is maximized; see Fig. 1(b). Then we constructed the matrix C and the vector B in 
(|4.3[) [see p.l4[[3~l~5)) ]. and neglecting the probabilities of having more than three photons inside the cavity, we assumed 
that we are given the following six frequencies v±i(m) (to = 1, 2, 3) normalized according to Y] T Y] ±j */ q (to) = 1. 
For simplicity we additionally assume that these frequencies are related as 

i^(l) = i^(2) = i/_i(2), i^(3) = i/_i(3). (4.6) 

For different values of i/i(l), z/i(2) and ^i(3) the numerical maximization of (|4.1|) over p a (m) under the constraint 
()4.3j) produced a result different from (14. 511 Q3 m =i S Q =±iP"( TO ) = -0- ^n example follows: for 

1^.(1) =0.05, i/i(2) = 0.25, i/i(3) =0.2 (4.7) 

the probabilities are: 

Pi(l) = 0.05148118, p-i(l) = 0.05087771, 
Pi(2) = 0.24811809, p_i(2) = 0.254403426, 

Pi(3) = 0.19158279, p_i(3) = 0.20353679. (4.8) 
Employing these probabilities in (|4.4[) and in (|3 . 9|) we get for the initial spin density matrix: 

p s = - [1 - (0.187183) 0-3 - (0.942992) a v + (0.275121) ctJ . (4.9) 

In this context we need to quantify the difference between the input frequencies v a (to) and the probabilities p a (m) 
which result from maximizing (|4. 1[) under the constraint (|4.3[) . In particular, this difference will quantify the relevance 
of the constraint (|4.3j) in maximizing (14. ip . A good measure of distance between two probability sets is provided by 

M 



6[u\\p] = 1 - E V^(m) Pa (m). (4.10) 



=±1 m=l 



This quantity is equal to its minimal value zero if (and only if) i/ a (m) — p a (m) (i.e., when the constraint (|4.3p holds 
automatically), and it is equal to its maximal value 1 for v a (m)p a (m) = 0. 

In Table I we calculated the distance 5[j/||p] between the frequencies and the corresponding probabilities. It is seen 
that in some cases this distance is just equal to zero, while for other cases it is rather small. 



5. CONCLUSION 



In this paper we describe a method for quantum state tomography. The usual way of solving this inverse problem 
of quantum mechanics is to make measurements of non-commuting quantities. Single apparatus tomography proceeds 
differently employing controlled interaction and measuring commuting observables. This is done via coupling the 
system of interest to an auxiliary system (assistant) that starts its evolution from a known state. The essence of the 
method is that the proper coupling is able to transfer the information on the initial state of the system to a commuting 
basis of observables for the composite system (system+assistant). 

It is important to implement the single-apparatus tomography for a situation with a physically transparent mea- 
surement base and with a realistic system-assistant interaction. Here we carried out this program for a two-level atom 
(system) interacting with a single mode of electromagnetic field (assistant). The atom- field interaction is give n b y the 
Jaynes-Cummings Hamiltonian, which has direct experimental realizations in quantum optics [l3l [iH IT?! Il8l. Il9l |20| . 
superconductivity p4| . semiconductor physics (25| . etc. As the measurement base we have taken presumably the 
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^i(l) = 0.05 


= 0.15 


= 0.25 


vi(l) = 0.30 


vi(2) = 0.05 


8 = 0.00989504 


8 = 0.0000494347 


8 = 1.1102 xl0~ 16 


8 = 0.00108428 


vi(2) = 0.15 


8 = 0.00318619 


5 = 


5 = 0.00140516 


8 = 0.0115233 


vi(2) = 0.25 


8 = 0.0000717018 


8 = 


8 = 0.0336704 




ui (2) = 0.30 


8 = 1.1102 xlO -16 


8 = 0.0000961022 







TABLE I: The distance <5[f ||p] (given by (|4.10[) 1 between the set of frequencies v±x{rri) and the set of corresponding probabilities 
p±i(m) obtained from maximizing (|4.1|l under the constraint (|4.3[1 . As in the main text, we assumed that the frequencies satisfy 
(|4.6[) . Thus the third frequency (not shown in the table) is got from /i a (3) = § — fJ>a(l) — /i Q (2). The numerical values for 
the involved parameters are: n = 2 (the initial average number of photons), g — 50 kHz (coupling constant), A = 100 kHz 
(detuning) and t = 300/^s (the interaction time). The matrix A4 and the vector B in this case are given by (|3.14l |3T5)) . Three 
places in the table are empty, because the corresponding frequencies are unphysical (their sum is larger than one due to the 
assumption (|4.6[) ). 

simplest set of observables related to the energies of the atom and field: population difference of the atoms a z and 
the number of photons a^a in the field. We have shown that one can determine the unknown initial state of the atom 
via post-interaction values of the average atomic population difference (172), the average number of photons (a^a) and 
the correlator of these quantities (& z a* a). The latter quantity does not need a separate measurement, since it can be 
recovered from the simultaneous measurement of the two basic observables a z and a^a. 

Since our scheme is based on measuring commuting observables, we can apply (more or less literally) the standard 
Maximum Likelihood setup for an approximate reconstruction of the unknown density matrix given the incomplete 
(noisy) measurement data. This is discussed in section 01 
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APPENDIX A: DERIVATION OF EQ. (f3TlO)) . 

Here we shall derive formulas for (cr z )t, (a^a)t and (a z a)a) t starting from (|3.1i 13 . 71 12. 9| . They are used in deriving 
(I3.10|) and they define the matrix M. in (|3.9p . One derives after straightforward algebraic steps 



-2ig(a x ) c ri 
n=0 
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ii in 



' n„ t ^ 
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where ((7j)o, i = x,y, z, are the unknown elements of the initial atom density matrix, {X + c.c.) stands for (X + X*), 
and where b\, c n and Q n are defined as 
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The average number of photons in the cavity, a) a, can be calculated in a similar way 

{a)a) t = 
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where 6 2 is defined as 



(A5) 
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The correlator of the two observables reads 

(a z <g> a t a) t = 
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where 63 is defined as 



& 3 = g — J^in + 1) [(2n + 3)c„+i - (2n + l)c„] 
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